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Abstract 

In this paper we study the dynamics in the general case for a Tavis Cummings atom in a non- 
uniform cavity. In addition to the dynamical Stark shift, the center-of-mass motion of the atom 
and the recoil effect are considered in both - the weak and the strong cavity atom coupling regimes. 
It is shown that the spatial motion of the atom inside the cavity in the resonant case leads to a 
transition between topologically different solutions. This effect is manifested by a singularity in the 
inter-level transition spectrum. In the non-resonant case, the spatial motion of the atom leads to a 
switching of the spin orientation. In both effects, the key factor is the relation between the values 
of the Stark shift and the cavity field coupling constant. We also investigate the entanglement of 
an atom in the cavity with the radiation field. It is shown that the entanglement between the atom 
and the field, usually quantified in terms of purity, decreases with increasing the Stark shift. 

PACS numbers: 73.23.-b, 78.67.-n ,72.15.Lh, 42.65.Re 
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I. INTRODUCTION 



Cavity quantum electrodynamics (CQED) is an active field of research focusing, on the 
quantum nature of the interaction of atoms with photons in high-finesse cavities UQ. Cur- 
rent issues of interest includes entanglement and quantum correlations js-S]. The standard 
model of CQED is a two level system in a quantum cavity. This case has served as an ex- 
ample for the realization of quantum programming protocols and for quantum teleportation 



10|. A controlled system of several atoms is also promising as a candidate for multiqubit 



entangled state. In the simplest situation of a single atom interacting with few photons in 
the regime of strong coupling, the coherent atom-photon interaction overwhelms incoherent 
dissipative processes. The system can then be described by the Jaynes Cummings ( JC) model 
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12| which captures the interaction of a single cavity mode with a frequency resonant 



with the transition frequency of the two lowest energy levels of the atom. The interactions 
between the atom and the radiation field should, however, involve not only the internal 
atomic transitions and field states but also should account for the center-of-mass motion of 
the atom and for recoil effects. Since the motion of the atom in a cavity and inter-level tran- 
sitions are connected to each other, the dynamics in a non-uniform cavity becomes rather 
complicated. For particular values of parameters, the corresponding classical system mani- 
fests regular or chaotic behavior: Different types of motion are found, including Levy flights 



and chaotic walking of an atom in a cavity 



13|. Consequences of the non-uniform cavity 



are also decoherence effects, a decay of entanglement and of the teleportation fldelity |14| . 
The main goal of the present work is the study of the dynamics of sm(1, 1) Tavis-Cummings 
(TC) system in a non-uniform cavity including the dynamical Stark shift (DSS) jisl . [l6| . 
This model describes two-photon transitions between the ground and the excited state via an 
intermediate state. The intermediate state can be eliminated from the equations of motion 



on the cost of introducing a dynamical Stark shift [15|. In particular, we study the impact 
of DSS on the dynamics of TC atom in non-uniform cavities. We identify different types of 
dynamics and possible mechanisms of switching between them. In the first part we use a 
semiclassical approach, which is a natural approximation for large mean photon numbers in 
the cavity. We investigate the problem in both the strong and the weak atom-cavity coupling 
regimes. In the second part we go beyond the semiclassical approximation and evaluate the 
infiuence of the cavity being non-uniform and also of DSS on the degree of entanglement 
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between the atom and the radiation field. 



II. MODEL 

The Hamiltonian of a single TC atom placed in an ideal cavity reads: 

H = —+ (a;o + C«^a)s^ + a;/a+a + 5((x)(s+a^ + s"(a+)^). (1) 

Here g{x) is the coupling constant between the atom and the radiation field, a, are the 
photon annihilation and creation operators, ( is the strength of the Stark shift, leading to 
the intensity-dependent transition frequency. The atomic two-level systems are described 
by the spin operators = , = ±i(T^), where (Jx,y,z are Pauh operators. If 

the initial kinetic energy of the system is small |^ <^ ^\ can neglect the atomic 



13| . Consequently, the operator 



2m ^ ""Y 2soV' 

motion and consider the standard TC model for uniform cavity g{x) ~ VLq. Here d is the 
atomic dipole moment, V is the volume of the cavity, Eq is the electric constant and n is the 
number of photons in the cavity. In the opposite limit the coupling constant g{x) depends 
on the position of the atom inside the cavity. Therefore, the motion of the atom has a 
strong impact on the inter level transitions leading to a complex nonlinear dynamics. We 
show below that this nonlinearity and the complexity can be utilized for a spin orientation 
control. We assume a standard form for the dependence gix) = cos kfX, where kj is the 
wave number in a lossless cavity of the Fabry-Perot type 
of the coordinate x together with the field and the spin operators form a complete set of 
observables. In what follows we will derive the Heisenberg equations for each observable 
and investigate the model ([1]) in the semi-classical approximation. We recall the standard 
commutation relations [s"*"; s~] = 2s^, [s^; s^] = ±s^, [a; a^] = 1 and the expressions for 
the operators in the interaction representation d{t) = de~^'^f^, d'^(t) = d'^e'^'^f^. From the 

Heisenberg equations of motion ^ = ^ + ^ [H; and eq. we infer 

dx 

— = ap, 
dr 

dp 

— = —smx-u, 
dr 

= -{5 + - + As.)v, (2) 

^ = (6 + — + As^u+ IGcosx ■ Sz ■ (N - + 2sl), 
dr ^ 2 ' 

dsz 

—— = — COSX ■ V. 
dr 



Here the semi-classical averaging procedure is used [13| x = kf{x), ay = {ay),ax = {ax),Sx 
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{sx), Sy = {sy), Sz = (sz), P = and the following notations are introduced a — — 



6 = -^Qp^, ~ ~ ^0^! u = A{axSx+aySy), v = A{axSy—aySx), a^+a^ — sl = N = const, 
sl. + s^+sl = const. We are interested in studying the system described by eq. ([2]) in different 
asymptotic limits. In particular, we focus on the dynamically induced control of the spin 
orientation and on switching. Our theoretical model can be realized easily in the experiment, 
e.g. by using Rydberg atoms and superconductive high-finesse cavities. For instance one 
may consider two-photon transitions between 44S'i/2 < — > 435*1/2 for rubidium atoms 85Rb, 
or 405*1/2 < — ^ 395i/2 for cesium 85 Cs. Such transitions proceed through the intermediate 
evel 39P3/2 or 43P3/2, and are of the type -two photons allowed and one photon forbidden- 
It! ]. Depending on the principle quantum number n the dynamical Stark shift is in the 
range < ( < lOOMHz [l^]. The cavity-atom coupling constant (which defines the time 
scale of our problem) is fio ~ 1 ~ lOMHz. Therefore, both limits > ( and Qq < ( which 
will be discussed below are realistic from the experimental point of view. 

III. ADIABATIC SOLUTIONS: THE RESONANT CASE 

In the semi-classical limit N ^ s^ and within the regime of strong coupling a ^ 1 eqs. 
([2]) tend to the form 

dx 

dp 
dr 

^ = -{g + Asz)v, (3) 
du 

— = (g + Asz)u+ 16 cos X ■ Sz ■ N, 
dr ' 

dsz 

— — = — cosx ■ V. 
dr 

From these relations it follows that the quantities 



W = -ucosx + ^{g + Aszf, dJid (4) 
u^ + v'^ + IQNsl = R^, 

are integrals of motion. Here g = 5 + ^ is the de-tuning parameter. 



— = — smx ■ ti, 
or 
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The set of equations ([3]) is characterized by two time scales. The small parameter a <^ 1 
sets the time scale of the slow center-of-mass motion of the atom; the fast time scale is 
determined by the inter-level transitions. The separation of time scales allows us to split 
the set of equations of motion into two parts 

— + au{T) sin X = 0. (5) 

and 

^ = {g + As,)u + 16cNs,, (6) 
^ = -cv. 

Due to the different time scales the slow parameter c = cos a; in eq.(l6]) hardly varies on a 
reasonable time scales. Since i; ~ a, a small change of the variable x(r) during the time 
interval t < 1/a can be neglected. The period of the inter-level transitions (see Eq.(12) 
below) is, however, shorter than T ^ 1/a. Consequently, the system performs a large 
number of inter-level transitions with c = cos(x(r)) ~ const. Utilizing the integrals of 
motion (jlj), one can derive a self-consistent equation for the spin operator Sz{t): 



dT 



-cv 




2Ac 



(9 



w 

c 



(7) 



Here in the energy of the system (4) we neglect the adiabatic part ^| ucos(x). Conse- 
quently in the resonant case, i.e. for (7 = we have: W = = + §, and eq. ([7]) 
can be rewritten in the form 



= -ci?oV(l-s2)(l 

A 



K 



■21 



K 



2cRq 



Ro = aVn. 



(8) 



Consequently, the spin dynamics is described by the following solution 



SAT 



cn{cRoT, k), k < 1, 
dn{cRoTK, i), K > 1. 

Here, cn{cRoT, k) and dn^cRrn, ^) are periodic Jacobi elliptic functions 



(9) 



. For the other 



variables in eq. (6) we obtain (upon straightforward but laborious calculations) the expres- 
sion 



K < 1, 



A sn'^{cRoT,hi) 

I ^sry{cRoTK, i), K > 1, 



(10) 



{Rosn{cRQT; K)dn{cRoT; k) k < 1, 
/ (11) 
^sn{cRoTK, K^sn^icRoTK, i) k > 1. 



From eq. ([9]) we conclude that, depending on the parameters of the problem, the dynamics of 
the level populations follow different solutions. These solutions are separated by the special 
value K = 1 of the bifurcation parameter k that signals the presence of topologically distinct 
solutions. In equation ([9]) the Jacobi elliptic functions are periodic in the argument with the 
period 

-^K(k), k<1, 
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where K{k) is the complete elliptic integral of the first kind 18|. If k — > 1, the period 



becomes infinite because K{k,) — > ln(4/vT^^ft?). The evolution in this special case is given 
by the non-oscillatory soliton solutions: 

Sz{r) = — ^ , 
cosh(citorj 

sinh^(ci?or) 

n^) = Ro — ,2. p ^ > (13) 

cosh [cRqT) 
u{t) = — — tanh^(ci?oT). 

The existence of a bifurcation parameter indicates that the solutions separated by it, 
have different topological properties. The phase trajectories corresponding to the solution 
dn{cRoTK, j:) are open and they describe a rotational motion, while trajectories correspond- 
ing to cn{cRoT, K,) are closed and they describe the oscillatory motion 19[. Due to the three 
independent variables in eq. (jH]) and the two integrals of motion (jl]), the system, described 
by eq. ([6]), is effectively one dimensional. Consequently, utilizing the integral of motion 
R^ = Rq = IQN , m(0) = 0, f (0) = 1, 3^(0) = 1 one can easily reduce the system (eqs. ([6])) 
to the effective one dimensional model 

Heff = '^ + ^^{l-slf + lRlsl. (14) 

It is straightforward to conclude that the solutions eqs.fl9|l-( fTT|) satisfy Hamilton's equations 
for taking the canonical pair of variables as (s^, v) 

^^ = -H. = -^, - = Hs. = ^. (15) 
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Due to the nonlinearity, the model expressed by eq. (HM shows a rich topological structure 
of phase-space trajectories. Topologically different types of phase trajectories are divided 
by bifurcation points which can be identified by evaluating the following Poincare index 

.^±/4ta„-.[-f]}. (16) 

Here 7 is a contour around the equilibrium point = 0,u = 0. From eqs. (fT^ . and (fT6|) 
we conclude for eq. (fT6i) in the linear approximation 

noRc 



cos'^(x) — ' — ^ 



J 



noRo 



(17) 



The change of sign in J from J = ltoJ = —1 marks the transition from the stable 
equilibrium point to the unstable equilibrium saddle point, i.e. from closed to open phase 
trajectories. The bifurcation point is defined by the simple relation x ~ arccos ^-^^g 
Obviously, the key issue is the relation between the position x{t) of the atom inside the 
cavity and the value of the dynamical shift of the frequency (. For the determination of the 
transition time on a quantitative level we need the exact solution of Eq. ([5]). For further 
analysis of eqs. ([5]), and ([6]) we note that both solutions, given by eq. iQ, contain the slow 
parameter k{t) = ( 2iioQois(a:(r)) )^^^- ^ initially ^^j^ < 1, cos(x(0)) ^ 1, and k{0) < 1, even 
in this case, due to the adiabatic motion of the atom inside the cavity, the condition k,{t) > 1 
can be realized as well. This means that the motion of the atom inside a non-uniform cavity 
leads to a tunneling of the system, formally through the presence of a separatrix. This leads 
to a qualitative change of the dynamics of Sz{t), i.e. eq. (|9]). The explicit solution for x{t) 
can be found from eq. ([5]). Taking eq. ( ITOj) into account and considering c = cosx as an 
adiabatic parameter on the time scale r ~ 1 /i?o from eq. (|5]) we infer 

^ + /(x) = /(x)cos(2i?or). (18) 

Here e = ct-^, f{x) = —etanx is the small parameter that controls the time scale for the 
slow motion of the atom inside the cavity and which is definitely larger than the time of 
the inter- level transitions \J^^ > (since > 1, a << 1). This means that during a 
slight change of the center of mass position from its initial value x{t = 0) = 0, the system 
performs multiple inter level transitions described by the first solution of given by eqs. (fTOj) . 



For solving eqs. ( !T8|) we make use of the existence of a slow and a fast time scale and look 
for a solution of the form 

x(r) = xi(r) + /iX2(r). (19) 

The time scale for the slow variable xi{t) is governed by 1/y/e. For the fast variable 2:2 (t) 
we have 2i?o, 2i?o ^ \/s and /i = With /(xi + yUX2) ~ /(a^i) + /^^2(§^)^^ "we infer from 
eq. (HI]) that 



dx 



+ 



f{Xi) +/iX2 



(9x 



cos(2i?o''")- 



(20) 



After averaging over the time interval l/2i?o we find 



Xl 



f{Xi) + ( /iX2 ( ^ ) cos(2i?o^ 



dx 



Xl 



fiX2 ~ —^X2 + f{xi) cos{2Rqt). 
Considering that (1x2 ~ fi{2Ro)'^X2 ~ 2_RoV^a^2 > P-X2i |^ ) eq. (122|) . The result is 



(21) 
(22) 



X2 r 



/(xi 



cos(2i?oT), 



.2/2ov^. 

Finally, in the linear approximation f{x) = —ex we conclude from eq. (I23ll that 



(23) 



x(r) = sinh 



1 + 



cos(2i?or) . 



(24) 



Using eqs. ([8]), and f l2^ we identify the time of bifurcation via 



|ln(arccos(^^)), k{rs) - I. 



(25) 



Eq. (l25l) connects the values of the DSS parameter C, and the birfurcation time r. This 
means that by observing the singularity in the spectrum of the Rabi oscillations of the 
inter-level transitions one can measure the frequency shift indirectly. Before proceeding 
with the non resonant case we summarize briefly the main results obtained for resonant 
case. The time scale of the spatial motion is described by the parameter a = Since 
X ~ a in the limit of a strong coupling between the atom and the cavity fio >> 7^? the 
motion of the atom inside of the cavity is adiabatic. Therefore, the system performs a large 
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FIG. 1: Color online. The time evolution of the atomic inversion, as descirbed by Eqs.dHl for 
the following the parameters: Rq = 20, a = 5 • 10~^, C/^o = 2. When the parameter k,{t) reaches 
the bifurcation values K^Tb) = ( 7^77-;; — ^-j—^ — 7^ ) ~ ''"ft ~ ( arccos (2^0^) ) ~ ^ the 



^ 2R0Q0 cos (a;(Ti,)) 

inversion oscillations manifest a singularity. 



number of inter-level transitions before the slow adiabatic coordinate changes significantly. 
If the resonant condition ( = 2{ujo — ujf) it met the inter-level transitions are described by 
topologically different solutions (cf. eq.(9)). The bifurcation parameter /c(a;(r)) depends on 
the center-of-mass coordinate of the atom x(r) (eqs.(24), (25)). We argue that the existence 
of topologically different solutions divided by a bifurcation point is the reason why the spatial 
motion of the atom leads to singularities in the spectrum of inter-level transitions. 



IV. ADIABATIC SOLUTIONS. NON RESONANT CASE 



In the non-resonant case fi' = ( § 



the form 



+ ojf — ojq^ ^ 0, the solution of the equation ([7]) has 



sAt) = 1 



1 



2W{cRt:^2;93)-A2' 

Here W{cRt; g2; gs) is the Weierstrass function [18| and the following notations are used: 



(26) 



g2 = 3Ai - 2Ai, gs = A1A2 - A\- -Aq, 



A2 = -^(l + /.^(2 ■ 2^^° 



A = -(2 + ^ 



(27) 



A, 



c 



2r2o-Ro cos(x) 
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cRr 



FIG. 2: Switching of the spin projection Sz{t) as induced by the motion of the atom inside the 



cavity. We have chosen 



Roc 



1. 



For K, « 1 i.e. C << 2fio-Ro, |^2| > |^i|, |^2| > |^o|, 92 = 3^4^ 5-3 = -A^ and the roots 
of the equation 

4x' - g2X -93 = 0, (28) 



satisfy the condition ei = —A2, 62 = 63 



2 ' 



^2 = 1(1 + 77^ I . Therefore we can use the 



foUowing representation of the Weierstrass function |l8| : 



W{t, g2, 93) = ei + ^ei cot^ ( V ^^1^ 




(29) 



As a resuh finally from eq. f l26p we obtain 

s.ir) = 1 



3ei + fei cot^ ^^J^icRt^ 



(30) 



From this expression we see that the motion of the atom leads to a switching of the spin 
orientation between the values sJt) = 1 and sJt) = 1 — tt-. The time interval for the 



3ei 



switching is set by the simple relation T 



An illustration of this switching is shown 



2cRo vBei ' 

by the simulations presented in Fig|2J This dynamically induced switching can be utilized 
for the manipulation of the spin orientation. A further important issue is the principle 
difference between the results obtained for the non-resonant case {( ^ Uf — uq) from those 
results corresponding to the resonant case {( = 2{uf — ojq)). In contrast to the resonant 
case, for the non-resonant situation the spatial motion of the atom leads to a switching of 
the spin orientation. However, due to the absence of a bifurcation and topologically different 
solutions, the singularities in the spectrum of inter- level transition are absent. 
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V. DYNAMICS FOR SMALL DSS AND MINIMAL CHAOS 



If the detuning between the radiation field and the spin sphtting is larger than DDS, i.e. 



for 6 = -7^^ ^ the equations of motion reduce to 

(Px 

1)-—- + au(r) sinx = 



-6ih 



2)( 



du 
dr 



dr 



COS XV. 



Since a <^ 1, the adiabatic approximation is still valid and the eqs. 
integrable. Introducing u(r) = X{t), v{t) = Y{t), Z{t) = RqSz, Cq 
solution can be written in the compact matrix form 

Y{t) = M{t) 
\Zir)) 



(31) 
(32) 

is analytically 
—RqCosx the 



^ x{o) ^ 

F(0) 



(33) 



M{t) 



5^ 



cosQnT — 



sinfiTvT 



■sinfiATT cos^att — ^sinfijvT 



Sf(l-cosl]^r) g^sinfi^vr 



^ + cos Q NT 



J 

6"^ = 6"^ + RqCos"^ X. For the particular initial conditions u{0) 
Sz{0) = 1, the solution simplifies and reduces to the compact form: 



where Qj^ 



CI 



(34) 

v{Q) = 0, 



U[T] 



IQN5 
02 

"TV 



(1 — cos{VLnt)) , 



V[T] 



N 



52 

-^ + --^cosQMr. 

N N 



— — — uj^{l — cos^tyt) sinx = 0. 



(35) 



(36) 



Using ( l35i) . equation ( 1311) can be rewritten as 

(fx 

Here 00 = |^ a/ a\6\. Equation (l36l) corresponds effectively to the perturbed universal Hamil 
tonian 



H 
Ho 

V{t) 



Ho + Vir 
1 



X^ — UJ^ cosx, 



(37) 



[ cos(x + ^nt) + cos(x — ^tyt)] . 
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FIG. 3: Color online. The chaotic motion of the atom inside the cavity in the regime of a weak 
coupling. The numerical integration of eqs. ([2]) is performed for the following values of the param- 
eters 6 = -0.5, a = 0.5, N = 50, s = 1, x{0) = 0, p(0) = l,n(0) = 0, v{0) = 0, s,{0) = 1. As 
evident from the figure, the motion in the nonadiabatic case is chaotic and resembles a diffusion 
process. 



This model shows a behaviour known as minimal chaos, as discussed in details in Ref. 19 1. 
The width of the stochastic layer formed near to the separatrix due to the time dependent 
perturbation V{t) can be estimated via the following expression: 

+00 +00 ^ oo 

AE = J {Hq;V} = — J ( sin(x + fi/vrr) + sin(x — r27vT))(ir = — y x sm(x — ^Ij^r). 

— oo — oo oo 

(38) 

Using the separatrix solutions for the unperturbed part of Hamiltonian Hq (I57|) . 19 | 



Xg 

X a 



we obtain 



4 arctan exp [ ± a;(r — Tq 
cosh \u}{t — To)] ' 

47rf2^ exp(7rr2Ar/2a;) 



< 



(39) 



(40) 



w sinh (^) 

Only in the energy interval located near the separatrix \E — Es\ < 6 Eg, the motion is 
irregular. However, in the non-adiabatic case, i.e. if a = is not small anymore, the 
dynamics is irregular in the whole phase space, as illustrated in Figs. (El IH andE]) 
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FIG. 4: Color online. The dynamics of the atomic inversion as deduced from the numerical 
integration of eqs.([2]) for the following parameters values: 5 = —0.5, a = 0.5, N = 50, s = 1, 
x(0) = 0, p(0) = l,n(0) = 0, ^(0) = 0, 5^(0) = 1. Obviously the inter-level transitions are of a 
chaotic nature. 



xlO' 

7\ ^ 
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FIG. 5: The correlation function for the atomic inversion as delivered by the numerical integration 
of eqs.© for 6 = -0.5, a = 0.5, iV = 50, s = 1, x(0) = 0, p(0) = l,u(0) = 0, v{0) = 0, s^(0) = 1. 
The finite width of the correlation function signifies the dynamical stochasticity. 



To confirm the existence of chaos, we examine the width of Fourier transform of the 
correlation function Gs^{t') = (s2(r + r')|s^(r)), Gg^ = J drCs^^r) exp[iuJT] = y+SM'" 

— oo ° 

+ 00 

With (...) = lim ^ f {. . .)dt we mean the averaging with respect to time, while Tc being 
the correlation time. The result of numerical calculations is presented in FigJS] The finite 
width of the Fourier transform signifies the emergence of chaos. The numerical results 
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presented in Figsj3l|5] confirm that in the hmit of weak couphng a = ~ 1 the dynamics 
turns chaotic. In contrast, in the hmit of a strong couphng chaos appears only in a smaU 
area near to the separatrix, as quantified by eq. lHOl) . For nonzero DSS the dynamics is 
integrable in both cases: for the resonant case ujq = uj + ^ and for the non-resonant case 
Wo7^w/ + | (EH]). 

Summarizing this part of the work, we found that when the detuning between the ra- 
diation field and the spin splitting is larger than DSS Uf — Uq > ( two different types of 
dynamics are realized: 1) In the regime of a strong coupling between the atom and the 
cavity a = ^ 1 the dynamics is regular in the whole phase space, except for a narrow 
energy gap near to the classical separatrix \E — Es\ < SEg (cf. Eqs. (39), Eq.(40)). 2) In the 
regime of a weak coupling a ~ 1 almost the whole phase space shows chaotic dynamics. Due 
to the impact of the spatial motion on the spin dynamics, the spectrum of the inter-level 
transitions is chaotic as well (see Fig.4). 

If the mean photon number in the cavity is not large, then the semi-classical approx- 
imation is not justified and the problem should be considered quantum-mechanically. A 
relevant question in this case is the quantum correlation within the compound quantum 
systems, i.e. between the atom and the radiation field. The entanglement between the 



atom and the field is usually quantified in terms of purity |2ll . 1221 ]. As entanglement is a 
specific quantum form of correlation it exhibits a number of essential differences to classical 
correlations. These issues will be the subject of the next section. 



VI. ENTANGLEMENT BETWEEN THE FIELD AND THE ATOM: QUANTUM 
MECHANICAL CONSIDERATION 

As was mentioned above, if the mean photon number inside the cavity is not a large 
number, then a semi-classical treatment is not viable and one should resort to quantum 



mechanics 
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22l |. Having said that, the translational degrees of freedom can still be 
treated classically; the inter-level transitions are described quantum mechanically, however. 
Thus, the center-of-mass coordinate of the atom x{t) may still be described by the classical 
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Hamilton's equations of motion 

X = ap, p=—{u)smx. (41) 

The quantum mechanical average should be performed then, i.e. 

(^) = (^(e^W + e-^^V(a+)^)y (42) 

The average is taken with respect to the wave functions that solve for the Schrodinger 
equation 

^^^ = H^nt\m)■ (43) 

Here 

Hint = (a+as, + g{x) {e'^'s+a^ + e-'^'s' {a+Y) , 6 = uJo-2ujf. (44) 

Following the standard procedure of Ref . [l^ , we seek solutions of the equation (143!) in 
the form 



\m) = E 



Ce,n{t)\e,n) + Cg,n+2{t)\g,n + 2) 



(45) 



Substituting eq. (1451) into eq. (143|) we find a set of equations for determining the expansion 
coefficients 

zCe,n = C^^^e,„, + ^7(x)v/(^rTI)(^r+2)e-*'^*C'3,„H_2, 

tCg,n+2 = -C^^C'<,,n+2 + 9ix)^/{n + l)(n + 2)e*^*Ce,n. (46) 

In order to derive analytical results for the degree of entanglement we will study the eqs. 
(H6|) in two different limits: For a strong atom-cavity coupling a = <^ 1, the quantity 
x{t) is an adiabatic variable on the time scale set by the inverse Rabi frequencies. In this 
case system, as described by eqs. ( l46l) . is exactly solvable and the solutions are given by the 
following expressions 



^ j (^-o* r ib{n) . /X rn\ , f r\ ^\ «(5 + C('^ + l))sin(A„t) ^ ■ 
Ce,„(t) = 2 <{ _ ^_^sm(A„t)Cg,„+2(0) + ( cos(A„t) - ^ ^ ^ )C'e,n(0) [>, 



C„„«(*) = - !M.i„(V)C,„(0) + (coB(A„i) + ''' + '^'"+^'""''"'> )c„„,.(0)}. 

6(n) ^ .^^^(ra^FliM. (47, 

Using these solutions and after tracing out the field's degrees of freedom one can introduce 
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FIG. 6: Color online. The purity P{t) as a function of time t and of the center-of-mass position of 
the atom, as dictated by eq. Eq.(51). The following parameters have been chosen Qq = 1 ,n = 1, 

C = i- 

the reduced density matrix for the atomic subsystem as 



Pn P?2 
P21 P22 



(48) 



where 



P12 — ^ Ce,n{t)C*g^n 



+21 



P2I ~ "y^j ^q.n+2(^)^e,n, 



Pit 



Cg,n+2(t) 



Ce,n(0) = CeWn, Cg,n+2(0) = = 1, = 0. 



(49) 



Here = is the distribution function of the field coherent states [12|. Substituting 

eq.( l471) in eq.( l49|) we conclude the following explicit form for the elements of the reduced 
density matrix ( HHj) (cjq = 2uf + (,n ':$> 1): 
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= i( 1 + e-«(i-cos(2U)) . cos(nsin(2tA) ) + 



2 , 
2 4A2 



-n(l-cos(2tA)) 



cos(?7, sin(2tA) 



= _!^e-*^*e-^(i-^°^(2tA))gi^^-gi^^2tA)+ (50) 

+ ie-C* . 4# ( 1 _ e-*^(i--(2tA)) . cos(nsin(2tA) 
2 2A^ \ 

= !^e*ae-^(i-'=°^(2*"))sin(nsin(2tA) + 
2A 

+ le<* . ( 1 _ e-..(cos(2a)-i) . eos(nsm(2tA) 

2 2A^ y 

P22 = 1 - e-^(i-™^(2tA)) . cos(nsin(2tA) 



2 A2 

i2 _ / S \ ,2/ 



The interaction between the atom and the radiation field is described by the nonseparable 
wave function (45), i.e. by an entangled state 12[. The entanglement between the atom 
and the field has a different meaning from the usual definition of the entanglement between 
the atomic states. The field is the subsystem with a large number of degrees of freedom and 
is usually prepared in a coherent state. Therefore, the state of the field is not infiuenced 
by the atom-field coupling interaction. For quantifying the entanglement between the atom 
and the field we should average and trace out the field's states on the cost of a partial loss 
of coherence. To make this point clear, let us consider the simplest protocol of a quantum 



measurement 



12| . The probability W(t, \ipait)), \ip jieidit))) that both subsystems are in the 



particular states \4'ait)) and lipjieidit)) is defined via the following relation: 

W{t, |^„(t)), l^f^^t))) = \{i;f,,i, I iija I m)\', (51) 



n j=e,g 

Here 

1pen{t) = Cenit), i^gn(t) = Cg,„+2(t). (53) 

As simple quantum measurement one may perform on the system is that, one observes the 
atomic state for an arbitrary field state. As we already have mentioned above, the field is 
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prepared in a coherent state I4]. The entanglement between the atom and the field (while 



statistically averaging over the field states) is quantified in terms of the purity P [lO| 



P{t, I ^a)) = E I E i^Jn{t){^a \ J,n) \' . (54) 
n=0 j=e,g 

Comparing eq. (54) with eq. (51) we conclude on a partial loss of coherence, since some 
interference terms are omitted in eq.(54). Nevertheless, the atomic and the field subsystems 
are still entangled and the density operator of the total system cannot be written as a direct 
product of density operators each corresponding to the atom and the field subsystems Using 
eqs. (49, 53 and 54) the purity of the quantum state is expressible in terms of the reduced 
density matrix in the standard form 10|, meaning that 

P{t, I = Pit) = Tr(^{p^Y^ = {p1,Y + (p^^)' + 2p?2P2i = 

= 1 - ^^{1 - exp [ - 2n{l - cos(2tA))] }, (55) 
A 

< P{t) < 1. 



From eq.(55) we readily deduce that the degree of coherence degrades with the decrease of 
DSS. The purity as a function of time and of the adiabatic coordinate is plotted in Fig.6. The 
purity is distributed non-uniformly, since the Rabi frequencies are coordinate dependent (H7I) . 
In the semi-classical limit, when the mean photon number in the cavity is large n ^ 1, the 

expression given by eq. (55) simplifies and for the time- and the coordinate-averaged purity 

(2^72 \ 
1 — ^2^4Q-2 1 • We see that the ratio of the cavity-atom coupling constant 

and the DSS fio/C is still the determining factor in our problem. From the expression above, 

the maximum of the purity is Pmax = 1 for the case Qq/C 1 and Pmin = | for ^o/C ^ 1) 

respectively. In the opposite case corresponding to the weak atom-cavity coupling with 

a = ~ 1, the analytical solution of the system (l46i) can be found in the special case of 

a resonant driving 6 = cjq — 2u!f = 

Ce,n = exp - (c^Q[nu{t)] + C2Q*[nu{t)] 

Cg,n+2 = exp ^!(!i±lM^ (^CQi^uit)] - C2Q*[nu{t)]^ , (56) 



where 



^ Ce,n(0) + g,,„+2(0) ^ Ce,.(0) - Cg,n+2i0) 

= ^ , ^2 = ^ , (57) 
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FIG. 7: Color online. The purity P{t) as a function of time t, as described by Eqs.(55, 57, 58). 
The following parameter are chosen: a = 10~^ n = 1, = 0.2, Qq = 0.1. 



and 



From eqs. (!56l [57|) we infer for the reduced density matrix and its purity the following forms 



Q[nu{t)] = exp in / u{t')dt' , Q ^[uj] = Q*[u], u{t) = cos{x{t)). (58) 



P{t) = Tr(^{p^Y^ = {pt,Y + (p^^)' + 2p?2P2i 
< P{t) < 1. 



(59) 



The explicit forms of the density matrix elements are 



Pl2 



P21 



P22 



n 

E 

n 

E 

n 



+ + C1C2 (g*' [na;(t)] + [na;(t)] ) 

- d + C1C2 (g*' [nu;(t)] - g' [na;(t)] ) 

- d + C1C2 (g' [nu{t)] - g*2 [r2a;(t)] ) 
d + Cl - C1C2 (g*' [na;(t)] + g' [na;(t)] ) 



(60) 



Considering the initial conditions (157|) . and tracing over the field states in eq. (!60!) . we 
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FIG. 8: The purity P{t) as a function of the time t. The results follow from Eqs.(55, 57, 61). In 
the calculations the following parameters are chosen a = 0.5 n = 1, = 0.2, = 0.1, oq = 1. 

obtain for matrix elements 

= lg-Ci| gnQ*[2u){t)]exp{-itC) _ gnQ[2a;{t)]cxp(-itC) jg-n^ 



4 
1 

4^ 



^ _gCU gnQ[2u;W]cxp(iiC) _ gnQ*[2a;{t)]cxp(itC) 



The time dependence of the matrix elements of the reduced density matrix, as given by eqs. 
(ISOj |6T|) . is governed by the exponential factors ( 158|) . We evaluate (5[2a;(t)] for two different 
regimes. In the regime of a regular motion, using solution ( l24l) we simply have 

Q\2uj{t)] ^exphiJ^lCiie'M) -Ci{l)]\, (62) 
L y «(, J 



where Ci{. . .) is the cosine integral function [18|]. In this case the purity exhibits fast 
oscillations, as evidenced by Fig. 7. From Fig. 7 we conclude that after a lap of time t = 
tc ~ 50, the oscillation frequency decreases. The reason for this behavior can be traced back 
to the time dependency of the parameters given by eq. (58) (i.e., Q[2uj{t)]): For t > t^, 
Q\2u(t)'j ~ exp [ — 2zy'^Cz(l)] = const and the diagonal elements of the density matrix 



20 



are constant in time 



Pii = 2 + 2 ^^P^ 



-n 



cos (nexp [2^ / —^Ci(l)] ) . 



(63) 



P22 = ^ - ^ exp(-n) cos {n exp [2 J ^Ci{l)] ) . 



(64) 



Hence, (pn) + (P22) — const and the temporal dependence of the purity is governed by 
the off-diagonal term 2p"2P2i- ^^e semi-classical limit n ^ 1 the purity turns constant 
P = 1/2 and is independent of the values of the frequency shift. In the regime when the 
motion is chaotic, we use a more advanced technique for the evaluation of Q[u!(t)]. Namely, 
because of the random character of the atomic motion, the exponent, given by eq. 0581) . 
should be considered as a functional of the random function uj{t). Therefore, we should carry 
out a statistical average with respect to all possible realizations of the random parameters. 
The mean value of the functional {Q[u{t)]) can be calculated by evaluating the following 
integral 



{QHt)]) = exp 



t 

i J uj{t')dt' 




lim 



/n 
duN ■ ■ ■ dui exp UkAtk Pn 





fc=i 



N 



1. 



(65) 



Here, the multidimensional normal distribution function is given via the expression 



Pn[uj] = (27r) 



d\\ . . . dX]\f exp - i ^ AfcWfc] exp - - ^ Ckk'^k^ 



1 



k,k' 



(66) 



where are distribution parameters and Ckk' is the covariation matrix 23[. Substituting 
eq. (166 p in eq. (16 5 p and performing the integration we obtain 



(g[a;(t)]) ^exp 




Erf 



(67) 



where Erf[. . .] is the error function, and 1/ao is the width of the correlation function 
c(r) = {uj(t + T)u{t)). Taking eqs. (161116 7p into account one can evaluate the time depen- 
dence of the purity in both the regular and the chaotic cases. The result of the numerical 
calculations is presented in Fig.dH]). As we see, after a fast decay, the purity stabilizes around 
the value P(t) ~ 0.5. The corresponding decay rate is determined by the correlation width 
l/ao of the random parameter x{t). 
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VII. CONCLUSIONS 



In the present paper, our aim was to analyze the dynamics of the TC atom in a non- 
uniform cavity, in order to uncover the consequences of the DSS for the weak and the strong 
cavity- field couplings. We found that in the regime of a strong coupling a = <^ 1, the 
motion of the atom inside the cavity is adiabatic. For the case when an exact resonance 
is reached ( = 2{ujo — Uf), then depending on the values of the center of mass coordinate 
of the atom x(t), the level populations are described by topologically distinct solutions 
separated by a bifurcation parameter. The motion of the atom inside the non-uniform 
cavity leads to a tunneling of the system through the separatrix and as a consequence to 
singularities in the inter-level transition spectrum Sz{t). The bifurcation point is identified 
by a simple analytical expression, given by eq.(25). A key factor is the relationship between 
DSS and the atom-cavity coupling constant C/^o- Therefore, by observing the singularity 
in the spectrum of the Rabi oscillations of the inter- level transitions Sz{t) one may infer the 
frequency shift in an indirect way. Far from the resonance ( » '2{u)q — uj), the spatial 
motion of the atom inside the cavity leads to a switching of the spin projection between 
the values Sz = I and Sz = 1/2. The period of switching depends on the values of DSS: 



T = 2fR..w.xi4-oV 9 = S: + In the regime of a weak coupling a = ~ 1, the 



motion of the atom inside the cavity becomes chaotic (see Fig. 3-5). We find that, in the 



adiabatic case, the ratio of the cavity-atom coupling constant and the DSS f2o/C ^-^'s still the 
determining factors. Namely, in the semi-classical limit, when the mean photon number in 



the cavity is large n >> 1, the time and the coordinate-averaged purity is P = ^1 — ^2_^_4Qi 
Therefore, the maximum of the purity Pmax = 1 is achieved for the case fig/C << 1 and 
Pmin = I corresponds to the case f2o/C >> 1- 
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